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A criterion for strong coupling between two quantum emitters and a single resonant light mode 
in a cavity is presented. The criterion takes into account the escape of cavity photons and the 
spontaneous emission of the emitters, which are modeled as two level systems. By using such 

criterion, the dissipative Tavis-Cummings ladder of states is constructed, and it is shown that the 
inclusion of one more emitter with respect to the Jaynes-Cummings (single emitter) case increases 
the effective parameter region in which n*'' order Rabi splitting is observed. 
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I. INTRODUCTION 

The study of light-matter interaction is one of the 
most fertile research areas in physics [1]. Under pre- 
cisely controlled conditions, matter and light can exhibit 
very interesting phenomenology. One such phenomena is 
the so-called Strong Coupling (SC) regime which is at- 
tained in the realm of cavity Quantum Electrodynamics 
[2]. This regime is achieved by isolating the matter and 
light in such a way that a cavity photon can interact 
several times with the matter forming an atom-photon 
"molecule" [3] . Experimentally, the interaction occurs in- 
side a cavity that modifies the spectral density of modes 
that the atom sees inhibiting spontaneous emission [4] 
and at the same time confining the light [5]. The forma- 
tion of this "dressed" light-matter state is experimentally 
verified in the photoluminescence spectrum of the system 
[6, 7]. As the bare matter and light frequencies become 
resonant, two peaks are observed [8, 9]. This is a manifes- 
tation of the formation of dressed states whose frequen- 
cies differ under ideal conditions by twice the radiation- 
matter coupling parameter, the so-called Rabi frequency 
[10]. Although this phenomenon is quantum mechani- 
cal, the Rabi splitting between the two modes and their 
widths can be obtained using a simple classical model of 
two damped harmonic oscillators [11]. The criterion for 
observing the Rabi splitting between one emitter and one 
mode is given by[ll-13]: 

5>|7-|, (1) 

where g is the light-matter coupling constant, 7± = 
(7a±7CT)/4, 7a is the decay rate of the photons of the cav- 
ity, which is inversely proportional to the quality factor 
of the cavity and 70- is the inverse of the lifetime of the 
emitter which is related to the Einstein A coefficient [10]. 
The condition given by (1) implies that the Rabi split- 
ting observed will be given by i? = \ — 7?. and that 
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the bare mode (light and matter) populations will oscil- 
late at the same modified Rabi frequency [12]. Although 
the Rabi spliting between the modes can be explained 
rising a classical model, the "quantumness" of the so- 
called linear Strong Coupling can be observed by other 
means, such as the counting statistics of the emitted pho- 
tons [11, 14, 15]. To evidence quantum behavior in the 
photoluminescence spectrum not affordable by classical 
oscillators, one has to consider the case in which more 
than one photon are bound to the matter. In this case, 
the emission frequencies exhibit anharmonicity [16] and 
thus, when the system cmiits a photon <:liffc;rc;nt emission 
peaks will be seen depending on the number of photons 
bound [17, 18]. The condition for observing Rabi split- 
ting between n photons and one emitter is now given by 
[13, 19]: 

V«5>|7-|, (2) 

and the Rabi splitting when there are n photons will 
be given by the Jaynes-Cummings (JC) modified Rabi 
frequency: 

Rn = ^j{V^g?~il. (3) 

This non linearity has been recently observed for a sin- 
gle emitter being an atom [20] or a superconducting 
qubit[21, 22] but remains to be evidenced in systems 
undergoing strong decoherence such as quantum dots 
in semiconductor microcavities. Much work has been 
devoted to the study of the emission spectrum of one 
quantum emitter coupled to a single mode and how de- 
coherence processes such as spontaneous emission, and 
finite life time of the photon, and different types of in- 
coherent pumping or finite temperature effects affect it 
[12, 17, 19, 23 26]. The case in which two atoms strongly 
interact with a single cavity mode has also been studied. 
The evolution of the closed systems was systematically 
investigated in [27-29]. The dynamics of the system in 
the linear regime was also investigated in [30, 31] and 
the evolution of the correlation functions as the system 
goes from one to several emitters was recently discussed 
in [32]. Finally, experimental results in which the linear 
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regime was explored for the case of two resonant excitons 
in separated quantum dots and whose Rabi constants are 
quite similar was presented in [33]. Nevertheless, no sys- 
tematic study has been presented of the dynamics in the 
non linear regime including decoherence. This is a rele- 
vant problem in quantum dots in microcavitites in which 
the effects of the environment on the quantum dot are 
quite strong and more importantly the quality factor of 
the cavities places the system in an intermediate regime 
[14, 19]. Understanding the dynamics of photons and 
qubits in a regime where the losses are comparable to 
the light-matter coupling, is relevant in several areas of 
quantum information processing in which the light mode 
acts as an information bus between processor qubits [34- 
39] and in general, as a test bed in which multipartite 
entanglement can be studied [40-43]. Finally, it is inter- 
esting to note that by preparing the initial state of the 
two emitters in their excited states one has direct access 
to the non linear regime even in the case where the cavity 
is empty and this might be easier to implement than the 
preparation of 1 or 2 photon states inside the cavity. In 
this work then, a criterion for observing strong coupling 
of two quantum emitters and a light mode in the non- 
linear regime in which anharmonicitics are apparent is 
presented. Although the results derived here are limited 
to identical quantum emitters where there is no incoher- 
ent pumping or finite temperature effects, they serve as 
the starting point for studying these extensions. 
The paper has been organized as follows: in section II, 
the Hamiltonian dynamics is studied, the eigenenergies 
and dressed states are written in terms of the coupling 
constant and the natural frequencies of the emitters and 
the mode; in section III, the Master equation that is 
used to model decoherence in the system is presented 
and the quantum regression theorem is introduced as the 
tool used to obtain the dynamics of the first order cor- 
relation functions of the subsystems; in section IV, by 
introducing the complex eigenenergies of the system, a 
criterion for SC is derived in terms of the emission rates 
of the system and the light matter coupling constant; fi- 
nally, some conclusions and general remarks are given in 
section V. 



the quantum emitters. The last term accounts for a 
dipole interaction between each quantum emitter and the 
light mode in the Rotating Wave Approximation (RWA) . 
The eigenstates of the Hamiltonian for the n^^ excitation 
manifold A„ in resonance (A = 0) are given by[42]: 

w(2/3) ^ _^ c,V4n - 2 

|a;(^)) = |n-l,5), 

where |n,j) = \n) \ j) in which the first ket is a 
Fock state of the field and the second ket is a Dicke 
state of the matter: |T_i) = |Gi) (Ga) , |To) = 
^(|Xi)|G2) + |Gi)|X2)),|ri) = \X,)\X,),\S) = 

^ {\Xi) IG2) - |Gi) 1X2)). The dressed state and en- 
ergy for the lowest excitation manifold is simply |0,T_i) 
with zero energy. For n = 1 the dressed states and ener- 

I (2)> I (3)\ I (4)> . (2) (3) (4) 

gies are \col } ,\u)l ) , \u)l ) with energies uj\ \ uj{ ' 
given by (5). Whenever the losses are very small, the 
emission spectrum of the system will consist of transi- 
tions in which the system goes from a dressed state with 
n excitations to another with n — 1 excitations, and the 
positions of the peaks in the photoluminescence spectrum 
will be given by the difference between their energies, 

i>"-"-i=a;W-a;^2r. (6) 

Finally, note if the emitters where not identical then the 
eigenenergies for the n**^ excitation manifold would be 
the roots of a quartic polynomial whose roots could only 
be expressed in terms of Ferraris solution[45] which are 
much complicated than the ones presented here and are 
of very limited use. 



II. HAMILTONIAN EVOLUTION 

The system is modeled using the Tavis-Cummings 
(TC) Hamiltonian [44] in which each quantum emitter 
has two possible states, ground \G) and excited \X) and 
the light is represented by Bosonic annihilation and cre- 
ation operators (a,a^). The Hamiltonian, in units in 
which h= 1, is given by: 

2 

H = woa^a ^ {(wq - A) ajai + giaja + a^^i)} (4) 

where ct^ = \Gi) {Xi\. The first two terms of the above 
equation represent the energies of the photonic mode and 



III. DISSIPATIVE DYNAMICS 

To fully account for the escape of photons outside the 
cavity and the spontaneous emission of the emitters, the 
dissipative dynamics of the system must be studied. This 
can be done by writing a Lindblad master equation for 
the density operator of the system that accounts for co- 
herent emission of photons (with rate 7^) and sponta- 
neous emission of the emitters (with rate 70-). Such mas- 
ter equation is given by [31]: 

|p = zip, H] + ^/:,{p} + ^ ^ (7) 

j=i 
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FIG. 1. (Color online) Ladder of bare states at resonance 
(A = 0) for two 2-level systems coupled to a cavity mode 
for the first three excitation manifolds and for two excitation 
manifolds containing n and n + 1 excitations. The double 
headed solid green horizontal arrows depict the light-matter 
coupling g, dashed red diagonal arrows the spontaneous emis- 
sion process 7a- and blue dotted vertical arrows the coherent 
emission process 7a. It is interesting to note that, because 
the spontaneous emission process acts independently on each 
two-level system, it is the only process that couples the singlet 
state IS') to the dynamics. The spacing between each excita- 
tion manifold A„ is given by ojo. In the non resonant case the 
atomic states with one excitation |n, 0) , |n, S) will be shifted 
downwards by an amount A and the state with zero atomic 
excitations |n, —1) will be shifted by an amount 2A. 

where Co{p\ = 20/50^ _ otQp _ pO'^O. A depiction of 
the processes included in the above master equation in 
the ladder of bare states is given in figure 1. 

To obtain the emission spectrum of the field or the 
emitters, it is necessary to write the first order correla- 
tion: 



c;o(i,r)-(Oni + T)0(t)) 



(8) 



where O can be any of cti, (T2, a. Once GQit, r) is at hand, 
the emission spectrum is given by [46]: 

S^{uj,T)^2kx (9) 



"'0 "'0 



where k is the finite bandwidth of the spectrometer and 
T is the amount of time in which light has been collected. 
To obtain the delayed time r dynamics of the correlation 
functions, the Quantum Regression Theorem (QRT) is 
used [47]. It asserts that, given a set of operators Xj 
satisfying the single time t dynamics. 



d_ 
dt 



(10) 



then the two-time dynamics with an arbitrary operator 
Y is given by: 

^ {X,{t + T)Y{t)} = -zJ2Ljk{Xk{t + r)Y{t)), (11) 

k 



for any operator Y. It can be easily seen that a^,cr| 
and crj can be written as linear combination of the ba- 



sis set, i^^j.™,,; 



dr 



\l,j){m,i\ with \l,j) e A„_i and 
|m, I) £ A„_i and that because the TC hamiltonian pre- 
serves the number of excitations N — a^a + ^I'^i 
they will satisfy the premise of the QRT. The expected 
values of the operators {Aj -.^ ■(t)) can then be arranged 
in a vector, x(i) that will satisfy (10) which can be 
written as: ^x(f) = — iLx(i) where L is a square ma- 
trix. If the two time expected values {Aj j.^ -{t + T)w{t)) 
{w e {(7i, (72, a}) are also arranged in a vector v{t + T,t), 
then, they will satisfy the same differential equation with 
respect to r, g|v(t + T,t) = — iLv(t + T,t). 

By ordering the operators Aj^.^^- according to the ex- 
citation manifolds they connect, the matrix L takes a 
block upper triangular form whose eigenvalues define the 
widths (Tj.) and positions (i>fc) of the emission peaks by 

Xk = -iffc + Vk- 

Note that, to use the QRT, the t — initial conditions 
are required; these are given by {Xk{t + T)Y{t)) \r=o = 
{Xk{t)Y{t)) = tr(p(t)Xfcy). The t dynamics of the ini- 
tial conditions required for the QRT can be studied in a 
similar fashion to the r dynamics. To this end, one sets 
T = to obtain the dynamics of the populations of each 
subsystem, Golt.T = 0) = {d\t)d{t)) = tr(p(t)Ol'd). 
In this case, it is easily seen that the required operators 
are linear combinations of the set Fij.„i i = \l,j) {m,i\ 
with \l, j) ^\m,i) € A„ and that this set also makes a 
closed set of differential equations. The expectation val- 
ues of such operators can be organized in another vec- 
tor y{t) that satisfies a differential equation of the form 
^y(^) — ~*Dy(i) where the matrix D is block upper 
triangular. 



IV. STRONG COUPLING CRITERION AND 
THE COMPLEX EIGENENERGIES 

In this section, the eigenvalues of the regression matrix 
L and the population matrix D are obtained and based 
on their dependence on the system parameters a criterion 
for observing Rabi splitting between the different transi- 
tions is derived. It is easily seen that L = wqI + B where 
B depends on g,ja,J<7 and A but not ujq, and 1 is the 
identity matrix. The Rabi splitting will then be deter- 
mined by the properties of L. In a similar way one can 
note that the population matrix D is also independent 
of Wo, and that the eigenvalues of both matrices can be 
written in terms of the complex eigenenergies in straight- 
forward generalization of equation (6). To this end in this 
section the complex eigenenergies of the Liouvillian (7) 
are introduced. In this case the real part will contain 
the information of the resonances of the system and the 
imaginary part will give information about its width. As 
we shall see the eigenenergies will still contain the sym- 
metries of the original Hamiltonian, i.e. they will come 
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FIG. 2. (Color online) Positions (thick lines) and widths 
(shadowing) of the complex eigenenergies at resonance for the 
first two rungs of the system as a function of the inverse life- 
time of the photon 7a. The spontaneous decay rate of the 
system has been set to zero (70- = 0). Note that in experi- 
ments typically luq ^ 9 and thus the widths between different 
excitation manifolds do not overlap. Also note that for more 
than one excitation there are always two eigenenergies that 
are degenerate but that differ in their associated width which 
is illustrated here for the second rung of the ladder. 



in a triplet and singlet that correspond to those in equa- 
tion (5) but now will have a non zero imaginary part 
that will account for the dissipation. Also as we shall 
show as the dissipation grows large the real part of the 
eigenenergies will be given only by nwo — A ± A and any 
trace of the coupling g will disappear. Note that even in 
the dissipative case the complex eigenenergy of the low- 
est energy-state, the light matter vacuum |0,r_i), is still 
strictly zero and is not affected by the dissipation. This is 
merely a reflection of the fact that the environment that 
was traced to obtain the master equation (7) is a zero 
temperature reservoir and thus the system will tend in 
the long time to reach its lowest energy eigenstate. For 
the first excitation manifold the complex eigenenergies 



are simply given by: 

Jl/2) 



= + Wo — A ± 7^1 

= -iy + Wo - A. 



(12) 



where the first order complex Rabi frequency is given by: 

(13) 



7^l = 1 2g2_ 7 



.A 



These are precisely the same energies that are obtained 
in the linear regime [30, 31], by treating the operators 
(Ti as bosons, and naturally reduce to the purely real 
eigenenergies of the Hamiltonian as 7a and 70- go to zero. 
From demanding that the modified Rabi frequency TZi be 
real at zero detuning the linear strong coupling condition 
is derived: 



V2g > |7_ 



(14) 



Jl/2) 



At zero detuning in the SC regime the widths of €{ ' are 
equal which stems from the fact that in the SC regime 
the dressed state picture is to some extent valid and be- 
cause of this at A = the two dressed states are half 
matter-half light and thus decay with an average of the 
two decay rates. Also note that the decay rate or width 
associated with e\ is purely due to 70- which also follows 
from the fact that the purely matter state |0,To) is an 
eigenstate of H. 

For more than one excitation (n > 0) the complex 
eigenenergies of the problem are simply given: 



.(1,2,3) 



- nwo - A + -P^I'^'S) 

- Tiwo — A, 
where the widths r„ are given by: 



2 

2 



(15) 



r„ = {n - l)7a -I- 7^ 



(16) 



and, 



Vl^^ = 7^„ cos 



l/cos(7r/6), (17) 



is given in terms of the discriminant 
6^3(7- +«f) 



and the complex Rabi frequency: 



/.9 



(18) 



' (4n - 2).g2 - 4 7_ + i 



(19) 



The complex eigenenergies are plotted in Figure (2) as a 
function of the inverse lifetime of the photon 7^. As it 
was mentioned before the eigenenergies still retain some 
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of the symmetries of the Hamiltonian, i.e. they arc spht- 
ted between a triplet and a singlet. Because the singlet 
|n — 1,5) whose energy is e^^ does not couple coherently 
to the other states of the system it simply acquires a non- 
zero imaginary part that corresponds to (n— 1) times the 
decay of a single photon plus the decay rate of a single 
atom. 

Far more interesting is what happens to the triplet of 
states. First note that the discriminant Q„ is a function 
only of the quantity 7_ + iA/2. This simply points to 
the fact that as far as the splitting between the different 
energies in a given excitation manifold is concerned, the 
decay rates act only as an "imaginary" detuning. 
To define the transition between strong and weak cou- 
pling it is necessary to study how the Rabi splitting at 
zero detuning (A = 0) between different energies in a 
given excitation manifold changes and in particular when 
does it become zero. Naively, it might be thought that 
the necessary and sufficient condition to have complex 
roots with non-zero imaginary parts is that the modified 
Rabi frequency, TZn be real (and thus automatically Q„ 
will be as well); this condition would be given by: 



\/4n - 2g > 2\j_ 



(20) 



Nevertheless, the condition given by the above equation is 
overly restrictive. Even in the case in which TZ^ is purely 
imaginary, the splitting given by (17) can have non zero 
real parts. Note that, when TZn = i^iTZn) is purely imag- 
inary, then also Qn = iS(Q„) and thus, the argument of 
the cos~^ in (18) is purely real; nevertheless, the cos~^ 
is a real number if its argument is real and in absolute 
value less than or equal to one. Thus even when TZn is 
purely imaginary if |3(Q„)| > 1 there will be n*^ order 
Rabi splitting. In summary, the most general condition 
for having non zero Rabi splitting in the n"' excitation 
manifold is, assuming that the complex Rabi frequency at 
zero detuning is purely imaginary, given explicitly by: 



|3(Qn) 



6\/3|7-/ff| 



(27_/5)2-(4n-2) 



> 1. (21) 



Again, note that if the Rabi frequency TZn is real at zero 
detuning then automatically the system is in SC. Only 
for cases when it becomes imaginary the above criterion 
is needed. In figure 3 the contour of |9(Q„)| = 1 is 
plotted as a function of n and ^-/g. It is clearly seen 
that whenever the contour crosses an integer n the Rabi 
splitting of the n*^ rung becomes zero. 
In the limit ^-/g <C 1, one can expand (17) to second 
order in 'j-fg: 

, 2 



■ 27_ 
*2n-l 

-5f\/4n^ 



2 + i 



2n- 



_ , _ (16n(n-l)-H) / 
-1 ''"i'23/2(2n-l)V2 y g J 



(22) 
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which explicitly shows that as 7^ and 7^ go to zero the 
complex eigenenergies become purely real and reduce to 



the eigenenergies of the Hamiltonian (4). 
Going back to the dynamics of the first order correla- 
tion function Q^(t,T), the matrix that characterizes its 
dynamics, L will be block upper triangular. Each block 
will represent the emission of one photon after the sys- 
tem decays from one excitation manifold to the one that 
is immediately below it. The matrix then will have a 3 di- 
mensional block representing the transitions between the 
first excitation manifold and the vacuum, then a 12 di- 
mensional block representing the transitions between the 
second and first excitations (with 4 and 3 eigenenergies 
that give 12 possible transition frequencies) and finally 
16 dimensional blocks representing transitions between 
the 4 possible energies of two contiguous excitation man- 
ifolds. The eigenvalues corresponding to the m^^ block 
will be given by 



(23) 



where z* is the complex conjugate of z. The last equa- 
tion tells that to obtain the transition frequencies the 
real parts are subtracted (in the Hamiltonian case this is 
precisely what is done (6)) and the imaginary parts are 
added to obtain the width of the emission line. Finally, 
it is worth mentioning that as can be seen from figure (2) 
at A = a pair of energies become degenerate and thus 
fewer peaks will be seen. For transitions between two 
excitation manifolds other than the one involving the 
vacuum, this will imply only nine transitions will give 
photons with different energies, and thus at most nine 
emission peaks could be observed. This number is fur- 
ther reduced if the initial condition of the system is such 
that its density operator has projections only in the sym- 
metric subspace spanned by the triplet states |n, Tj). In 
this case the system will not be able to transit through 
the singlet state \n, S) and thus the contribution asso- 
ciated with the energy will not appear in equation 
(23). As for the dynamics of the occupation numbers 
which serve as the initial conditions for the QRT the ma- 
trix D that characterizes its dynamics will also be block 
diagonal but know the blocks will correspond to energy 
differences between the same excitation manifold: 



rn rn \ m ) 



(24) 



Note that for the excitation manifold that corresponds 
to the vacuum there is only one energy eo = and only 

one (5o = which simply accounts for the conservation of 
probability (tr(/5(i)) = 1) that the master equation (7) 
provides. 



V. CONCLUSIONS 

A criterion for observing Rabi splitting between two 
quantum emitters and a single cavity mode has been 

presented in terms of the Rabi splitting of the complex 
eigenenergies of the master equation (7). It has been 
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FIG. 3. (Color online) The black thick line is the contour 
of |S5(Qn)| = 1 at resonance A = , n is taken as a real 
number whose scale is in the right of the figure. The purple 
thin lines are the different Rabi splittings j5R(en — ^Qn)\ at 
resonance A = for the first 4 excitation manifolds (the scale 
of the different splittings is in the left of the figure). Note that 
whenever the contour of |9(Qn)| = 1 crosses an integer value 
of n as 7-/(7 is varied one of the Rabi splittings becomes zero. 
Finally, as expected from the purely Hamiltonian picture the 
vF^ order Rabi splitting reduces to \/4n — 2 as the dissipation 
effects become negligible 7_ j g — >• 0. 

shown that the criterion given by (21) is a robust char- 
acteristic of the dynamics given by (7) since it accounts 
both for the delayed dynamics, which is relevant for the 
calculation of the first order correlation function and also 
the dynamics of the density matrix populations which 
contains the information about the occupation numbers 
of the subsystems. The precise positions of the emission 
peaks were given in terms of the real parts of the com- 
plex eigenenergies and the rich multiplet structure of the 
photoluminescence spectrum was discussed. An intuitive 
picture of the possible types of multiplets was given in 
terms of whether each rung of the Tavis-Cummings lad- 
der exhibited Rabi splitting or not. Although the condi- 
tion given by (21) determines the presence of oscillatory 
frequencies in the dynamics of the first order correlation 
function that will lead to anti-crossings in the photolumi- 
nescence spectrum, this anti-crossings will not necessarily 



be easily resolved as can be seen in figure 2. The reason 
for this is that the broadening of the spectral lines, which 
are given by the imaginary parts of the eigenenergies (15), 
grows at least linearly in "nqa- On the other hand, the 
spacing between the lines grows approximately as \/ng 
as can be seen from (22). Because of this last obser- 
vation, it is more feasible to observe the anharmonicities 
described here by focusing on suppressing the emission of 
the cavity mode, i.e., in very high Q cavities. The results 
presented here allow to analyze the differences between 
the well known single emitter (JC) case in the dissipative 
regime and the case were two emitters are present (TC). 
It was already well known that the number of emission 
lines increases significantly; this is simply because each 
excitation manifold will have more states and thus more 
transitions can occur. With the results presented here, 
also the dissipative open system dynamics can be com- 
pared. For instance, in the JC model, all the dynamics 
is determined by whether the modified Rabi frequency 
(3) of the one emitter case is real or not, whereas in the 
case of two emitters, a more involved criterion is neces- 
sary (21). The system with two emitters is more robust 
against decoherence since the parameter region in which 
SC can be observed is bigger as compared to the area in 
which the two emitter Rabi frequency is real. It is also in- 
teresting to note that, in the JC case, 7_ only shrinks the 
Rabi splitting between the lines but does not modify the 
broadening of the spectral whereas, in the two emitter 
case, it does, as can be seen from (22). Thus, this work 
gives insights in the interplay between cooperativity and 
dissipation by analyzing the simplest case of cooperative 
effects in the interaction between light and matter under 
decoherence effects. 
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